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Mathematical  modelling  is  used  to  analyse  water  management  in  Alkaline  Direct  Methanol  Fuel  Cells 
(ADMFCs)  with  an  anion  exchange  membrane  as  electrolyte.  Cathodic  water  supply  is  identified  as  one  of 
the  main  challenges  and  investigated  at  different  operation  conditions.  Two  extreme  case  scenarios  are 
modelled  to  study  the  feasible  conditions  for  sufficient  water  supply.  Scenario  1  reveals  that  water  supply 
by  cathodic  inlet  is  insufficient  and,  thus,  water  transport  through  membrane  is  essential  for  ADMFC 
operation.  The  second  scenario  is  used  to  analyse  requirements  on  water  transport  through  the  mem¬ 
brane  for  different  operation  conditions.  These  requirements  are  influenced  by  current  density,  evapo¬ 
ration  rate,  methanol  cross-over  and  electro-osmotic  drag  of  water.  Simulations  indicate  that  water 
supply  is  mainly  challenging  for  high  current  densities  and  demands  on  high  water  diffusion  are 
intensified  by  water  drag.  Thus,  current  density  might  be  limited  by  water  transport  through  membrane. 
The  presented  results  help  to  identify  important  effects  and  processes  in  ADMFCs  with  a  polymer 
electrolyte  membrane  and  to  understand  these  processes.  Furthermore,  the  requirements  identified  by 
modelling  show  the  importance  of  considering  water  transport  through  membrane  besides  conductivity 
and  methanol  cross-over  especially  for  designing  new  membrane  materials. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Recently,  alkaline  fuel  cells  (AFCs)  attract  attention  by  reason  of 
their  ability  to  operate  with  non-precious  metals  as  catalyst.  Due  to 
higher  stability  and  activity  of  non-precious  metals  in  alkaline 
media  compared  to  acidic  media,  it  is  possible  to  use  metal  catalyst 
in  AFCs  other  than  platinum  (Pt),  which  is  commonly  used  as 
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catalyst  in  low  temperature  fuel  cells.  Replacing  Pt  catalyst  by  other 
metal  catalyst  can  reduce  cost  of  fuel  cells  and,  thus,  is  the  main 
motivation  to  improve  AFCs.  Hence,  research  regarding  AFC  mostly 
focusses  on  new  catalyst  materials  for  oxygen  reduction  and 
oxidation  of  various  fuels  such  as  hydrogen,  methanol,  or  other 
alcohols  in  alkaline  media.  A  review  about  new  catalyst  for  AFC  can 
be  found  in  Refs.  1,2  . 

Alkaline  fuel  cells  are  commonly  operated  using  a  liquid  elec¬ 
trolyte  such  as  potassium  hydroxide  (KOH)  or  sodium  hydroxide 
(NaOH)  which  reacts  with  dissolved  carbon  dioxide  (CO2)  according 
to  Eq.  (2)  forming  carbonate  ions  (CO3-): 

CO2  +  2OH  +H20  (1) 


912 


C.  Weinzierl,  U.  Krewer  /  Journal  of  Power  Sources  268  (2014)  911—921 


2K+ +C0f-->K2C03|  (2) 

On  one  hand,  this  carbonation  process  decreases  the  perfor¬ 
mance  of  the  fuel  cell  due  to  slower  kinetics  of  fuel  oxidation  at 
anode  [3]  and  lower  ionic  conductivity  of  the  electrolyte  [4  .  On 
the  other  hand  carbonates  in  supersaturated  solutions  can  form 
precipitating  salts  (Eq.  (2))  that  may  block  pores  of  gas  diffusion 
layer,  catalyst  layer  or  membrane.  In  order  to  avoid  the  disad¬ 
vantages  caused  by  carbonation  of  liquid  electrolyte,  anion  ex¬ 
change  membranes  were  invented  as  a  solid  alkaline  electrolyte 
[5].  It  has  been  reported  that  similar  performance  of  a  fuel  cell 
can  be  obtained  with  an  alkaline  anion  exchange  membrane 
(AAEM)  in  CO3  form  and  with  the  same  AAEM  in  OH-  form  [6], 
Still,  there  are  only  few  studies  about  AFCs  using  solely  a  solid 
electrolyte  membrane  instead  of  liquid  electrolyte  and  most  of 
them  analyse  cell  performance  for  new  membranes,  ionomers  or 
catalysts  [7-11  . 

It  has  already  been  noticed  that  water  transport  through  mem¬ 
brane  might  be  a  limiting  factor  [11].  Nevertheless,  only  little 
attention  is  given  to  process  engineering  issues  like  water  transport 
through  the  membrane  or  water  management  so  far.  Diffusion  co¬ 
efficient  and  electro-osmotic  drag  coefficient  of  water  through  an 
anion  exchange  membrane  were  determined  from  experimental 
results  by  Ref.  [12].  Likewise,  the  water  diffusion  coefficient  through 
an  anion  exchange  Morgan  ADP  membrane  was  estimated  from 
sorption  kinetics  by  Ref.  [13  .  A  model  based  prediction  of  the  water 
diffusion  coefficient  through  a  SnowPure  Excellion  1-200  anion  ex¬ 
change  membrane  is  given  in  Ref.  [14]  and  electro-osmotic  water 
drag  coefficients  in  alkaline  media  were  also  determined  [15  .  Water 
management  has  been  analysed  for  AFCs  with  liquid  electrolyte  by 
Refs.  16,17  .  These  two  studies  are  based  on  mathematical  model¬ 
ling  and  belong  to  the  few  modelling  studies  that  consider  water 
transport  or  management  in  AFCs  [16-21  .  Although  water  man¬ 
agement  in  AFCs  is  even  more  challenging  without  liquid  electrolyte 
due  to  the  lack  of  the  buffer  function  of  alkaline  solutions,  water 
management  has  only  been  studied  for  the  anode  of  anion  exchange 
membrane  fuel  cells  so  far  [20,21  .  These  two  studies  are  using 
stationary  20  and  dynamical  21  3D  modelling  to  analyse  flooding 
of  anode  including  the  influence  of  changes  in  various  operation 
conditions.  Some  modelling  works  about  AFC  contain  analysis  of  fuel 
cell  performance  [22-26].  The  latter  is  the  only  modelling  study 
regarding  alkaline  direct  methanol  fuel  cells  (ADMFCs).  Further¬ 
more,  another  two  model  based  studies  regarding  AFCs  without 
liquid  electrolyte  were  found  in  literature,  both  of  them  analysing 
the  conductivity  and  water  uptake  of  the  AAEM  [27,28  .  Hence,  there 
is  a  lack  of  knowledge  and  publications  in  literature  regarding 
modelling  and  simulation  of  anion  exchange  membrane  fuel  cells, 
particularly  related  to  water  management  at  cathode  and  in  the 
whole  fuel  cell.  It  is  still  unclear  in  which  way  and  intensity  water 
management  influences  and  limits  fuel  cell  performance  and  what 
can  be  done  to  improve  water  management  and  to  widen  the 
operation  range  of  the  fuel  cell.  The  present  study  fills  that  gap  by 
investigating  limitation  of  current  density  by  water  management. 
Requirements  on  operation  conditions  and  water  transport  are 
analysed  in  order  to  estimate  conditions  for  stable  operation.  To  our 
knowledge,  this  paper  presents  the  first  model  based  analysis  of  the 
water  management  in  an  ADMFC. 


2.  Functional  principle  of  an  ADMFC 

As  all  fuel  cells,  ADMFCs  convert  chemical  energy  stored  in  a  fuel 
directly  into  electric  energy  using  electrochemical  reactions.  In  the 
case  of  an  ADMFC,  methanol  is  used  as  fuel  which  reacts  with 


hydroxide  ions  (OH~)  in  the  electrochemical  oxidation  reaction  at 
the  anode  producing  water  and  carbon  dioxide  and  releasing 
electrons: 

CH3OH  +  6OH-— >C02T  +5H20  +  6e-  (3) 

At  the  cathode,  OH^-ions  are  produced  in  the  electrochemical 
reduction  reaction  of  oxygen: 

|o2  +  3H20  +  6e-^60H“  (4) 

Hence,  the  overall  reaction  is  the  same  as  in  acidic  direct 
methanol  fuel  cells: 

CH3OH  +  |o2^C02  +  2H20  (5) 

During  operation,  electrons  are  transported  from  anode  to 
cathode  via  an  electric  circuit  whereas  OH_-ions  are  conducted  by 
an  electrolyte  from  cathode  to  anode. 

Since  C02  is  permanently  produced  at  the  anode,  the  carbon¬ 
ation  process  is  more  intense  in  ADMFCs  than  in  alkaline  hydrogen 
fuel  cells.  Hence,  it  is  important  to  replace  the  liquid  electrolyte  by 
a  solid  electrolyte  membrane  in  ADMFCs.  A  schematic  of  an 
ADMFC  with  a  solid  electrolyte  membrane  is  shown  in  Fig.  1.  The 
anode  is  fed  with  a  methanol  water  solution  while  the  cathode  is 
flushed  with  air.  Although  the  reactants  are  separated  by  a 
membrane,  some  mass  transport  through  the  membrane  takes 
place  in  addition  to  the  ionic  transport.  Methanol  transported 
through  the  membrane  by  diffusion,  hereafter  named  as  methanol 
cross-over,  is  oxidised  at  the  cathode  in  a  chemical  reaction  similar 
to  the  overall  reaction  Eq.  (5)  producing  water.  This  loss  of 
methanol  decreases  the  efficiency  of  the  fuel  cell.  Furthermore,  the 
methanol  oxidation  at  the  cathode  causes  a  decrease  in  cell  po¬ 
tential  due  to  a  mixed  potential  [29].  Similar  to  methanol,  water  is 
also  diffusing  through  the  membrane  from  anode  to  cathode  but 
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Fig.  1.  Schematic  of  an  ADMFC  including  anode  and  cathode  fluid  channels,  gas 
diffusion  layers  (GDLs),  catalyst  layers  (CLs)  and  a  solid  electrolyte  membrane. 
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there  is  also  a  back  transport  by  electro-osmotic  water  drag  caused 
by  the  ionic  transport.  Hence,  water  is  consumed  at  the  gas 
flushed  cathode  due  to  the  reduction  reaction  Eq.  (4)  and  electro- 
osmotic  water  drag.  Thus,  water  has  to  be  supplied  to  the  cathode 
to  provide  sufficient  amount  of  water  for  the  electrochemical  re¬ 
action.  At  the  same  time,  water  is  added  to  the  liquid  flow  at  the 
anode  due  to  the  electrochemical  oxidation  reaction  Eq.  (3)  and 
the  water  drag. 

All  these  information  lead  to  the  conclusion  that  water 
management  in  ADMFCs  is  challenging  at  both  electrodes  of  the 
fuel  cell.  This  paper  focuses  on  water  management  to  enable 
stable  cathode  performance  by  analysing  one  big  challenge 
regarding  water  management  in  ADMFCs:  the  cathodic  water 
supply. 


3.  General  mathematical  model 

In  this  study,  mathematical  modelling  is  used  to  describe  and 
analyse  processes  and  water  management  in  an  ADMFC.  For  that 
reason,  the  following  general  model  is  derived  which  is  modified 
later  forming  different  scenarios.  The  two  scenarios  analysed  in  this 
paper  solely  describe  a  single  fuel  cell  which  is  sketched  in  Fig.  2. 
The  membrane  which  couples  anode  and  cathode  is  not  modelled 
in  detail,  but  it  is  implemented  as  a  semipermeable  wall  that  allows 
transport  of  certain  species.  In  order  to  focus  on  the  important 
processes,  the  model  is  kept  as  simple  as  possible  by  including  the 
following  assumptions. 

Pressure  and  temperature  at  anode  and  cathode  are  assumed  to 
be  constant  at  1  bar  and  323.15  K,  respectively.  Therefore,  it  is  not 
necessary  to  include  an  energy  balance.  Furthermore,  charge  bal¬ 
ances  are  not  considered  since  the  model  is  not  used  to  study  the 
performance  but  to  analyse  procedural  requirements.  Hence,  the 


differential  equations  of  this  model  arise  from  mass  balances  of  the 
components  j3  with  the  space  direction  k. 


dpg  d  (  .  \  * 

IT  “  ~dTk  \peVk  +JkA  + 


(6) 


The  model  is  used  to  analyse  overall  water  management  of  the 
whole  fuel  cell  in  order  to  investigate  conditions  that  are  at  least 
required  for  sufficient  water  supply  to  the  cathode.  The  minimum 
requirements  can  be  obtained  by  assuming  an  ideal  system  in 
which  transport  processes  within  the  chambers  e.g.  through  gas 
diffusion  layer  (GDL)  and  catalyst  layer  (CL)  are  disregarded.  Thus,  it 
is  assumed  that  no  local  gradients  appear  in  anode  and  cathode 
chamber  which  are  modelled  as  continuous  stirred  tank  reactors 
(CSTRs).  The  consequences  of  this  assumption  are  discussed  for 
each  scenario  separately.  Furthermore,  diffusion  through  mem¬ 
brane  is  described  by  Fields  law.  Methanol  (Me)  diffusing  through 
the  membrane  is  immediately  oxidised  at  the  cathodic  catalyst 
layer  as  well  as  liquid  water  is  immediately  consumed  or  evapo¬ 
rated  at  the  cathode  (c£liq  =  0).  Water  drag  takes  place  from 
cathode  to  anode  (n^)  while  it  is  assumed  that  no  methanol  is 
dragged  by  ions.  Thus,  the  mass  balances  for  all  components 
^e{Me,  H20,  Meg,  H2Og, C02, 02,  N2}  reduce  to: 


V* 


V ( 


HCA 

•A 

•A 

dt 

-  in 

out 

dcft 

•c 

•C 

dt 

-  ni3, in 

tt/3,  out 

•diff  A 

ne,  +  an 


i  Adiff  i  J 

+  rip  +  a, 


(7) 


(8) 


with  li 


(3,  out 


A  pi 
L(Srout 


(9) 


Fig.  2.  Sketch  of  the  general  model  of  an  ADMFC.  The  chambers  are  well  mixed  (no 
local  gradients)  and  the  membrane  allows  transport  of  certain  species. 
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Neglecting  mixing  volumes,  the  volume  flows  leaving  anode  and 
cathode  can  be  derived  from  the  overall  mass  balances  in  the 
chambers: 


cA  _  cA 

^out  —  Mn 
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with  n^e  =  ri{J}e  and  c*/t  being  the  concentration  of  the  pure 
substance  of  component  a  at  anode  or  7  at  cathode  and 
a£{Me, H20, Meg, H2Og, C02}  as  well  as  Te{H2Og, C02, 02, N2}. 
The  inlet  flow  rates  at  anode  and  cathode  are  defined  by  con¬ 
sumption  of  the  reactant  proportional  to  an  excess  ratio  A  which  is 
kept  constant  in  order  to  avoid  reactant  starvation.  Consequently, 
volume  flow  rates  change  with  current  density  (Eqs.  (A.20)  and 
(A.21)). 

Since  the  model  is  not  used  to  analyse  dynamics  or  system 
performance,  the  kinetics  of  reactions  are  disregarded  and  reaction 
rates  are  described  by  Faraday's  law.  It  is  further  assumed  that  gas 
and  liquid  at  anode  are  always  in  equilibrium.  For  the  calculation  of 
the  molar  fractions  of  water  in  the  gas  phase  yHi 0g,  the  laws  of 
Raoult  and  Dalton  are  used  as  well  as  the  ideal  gas  law  and  the 
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Wagner  equation  [30  .  The  corresponding  equations  are  declared  in 
A.  Thus,  the  resulting  sources  and  sinks  ap  at  anode  and  cathode  in 
Eqs.  (7)  and  (8),  which  are  caused  by  reaction,  evaporation  and 
electro-osmotic  drag,  can  be  expressed  by: 


AMi _ y^eS  AMI 

W  ”1-J/MeS-J/H2OS  W 


(13) 


Table  2 

Physical  constants  and  geometry  parameters  used  for  simulation  of  all  scenarios. 


Parameter 

Value 

Parameter 

Value 

F 

96,485  A  s  mol-1 

R 

8.31  J  mol-1  K1 

VA 

12.5  ml 

V0 

12.5  ml 

am 

25  cm2 

dM 

30  pm 

5AMi 
6  F 


^H,Og 


1  yMes  ^H,0” 


AMi  AMi 

~6F+K~T 


A  _  AMi 
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(17) 

(18) 

(19) 


ADMFCs.  In  this  paper,  conditions  for  sufficient  water  supply  to  the 
cathode  are  studied.  Therefore,  the  following  two  scenarios  are 
created  that  fulfil  cathodic  water  demand  in  various  ways.  In  order 
to  keep  the  models  as  simple  as  possible,  these  scenarios  only 
describe  a  single  fuel  cell  fed  with  a  methanol  solution  of  1  mol  1  1 
and  a  constant  excess  ratio  Aa.  Hence,  the  molar  flows  entering  the 
anode  are  time-independent,  nAin  =£f(t)  but  do  change  with  current 
density.  Both  scenarios  are  simulated  using  Matlab.  In  the  first 
scenario,  the  algebraic  equations  displayed  below  are  calculated 
straightforward.  In  order  to  be  more  flexible  in  extending  and 
modifying  the  model,  the  second  scenario  is  numerically  integrated 
for  various  current  densities  with  the  solver  odel5  s  of  Matlab  until 
steady  state  is  reached.  The  steady  state  results  are  shown  in  this 
study.  These  results  are  same  if  solver  ode45  of  Matlab  is  used  for 
the  numerical  integration.  However,  it  would  also  be  possible  to 
calculate  the  steady  state  concentrations  by  solving  a  complex 
system  of  algebraic  equations  which  leads  to  the  same  results  as  the 
numerical  integration. 

In  the  first  scenario,  water  is  supplied  by  humidifying  cathodic 
inlet  gas.  This  scenario  is  used  to  analyse  which  inlet  conditions 
lead  to  sufficient  water  supply.  In  the  second  scenario,  water  is 
supplied  by  mass  transport  through  membrane  to  analyse  the  re¬ 
quirements  on  water  diffusion  for  sufficient  water  supply. 
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Further  equations  that  complement  the  model  are  listed  in 

Appendix  A. 

The  parameters  and  initial  conditions  listed  in  Tables  1-3  are 
used  for  the  simulation  of  both  scenarios  unless  specified 
differently. 


4.  Modelling  and  analysis  of  different  scenarios  for  cathodic 
water  supply 


4.2.  Scenario  1  -  water  supply  by  cathodic  inlet 

This  scenario  is  used  to  identify  inlet  conditions  that  lead  to 
sufficient  water  supply.  Therefore,  a  mathematical  model  is  pre¬ 
sented  first,  followed  by  the  discussion  of  the  results. 


4.2.2.  Mathematical  modelling 

In  order  to  study  the  effect  of  water  supply  by  solely  humidi¬ 
fying  the  cathodic  inlet  gas,  molar  flows  through  the  membrane  are 
assumed  to  be  zero  (rij^o  =  rij^5  =  nj^e  =  0).  This  is  implemented 
by  setting  the  corresponding  coefficients  equal  to  zero:  k  =  0,  D^1  = 
0  for  P  e  {Me,H20}.  Hence,  anode  and  cathode  are  decoupled  and  it 
is  sufficient  to  consider  only  the  cathode  for  the  analysis  of  this 


The  general  mathematical  model  above  is  modified  to  describe 
different  scenarios  for  the  analysis  of  the  water  management  in 


Table  1 

Chemical  data  used  for  simulation  of  all  scenarios. 


Parameter 

Value 

Parameter 

Value 

^Me 

-8.54582 

M\ Me 

32.04  g  mol-1 

#Me 

0.67266 

Mh2o 

18.02  g  mol'1 

ClVIe 

-2.54743 

mco2 

44.01  g  mol-1 

-2.71874 

M02 

32.00  g  mol-1 

^H20 

-7.71374 

LMe 

24.66  mol  1_1 

Bh2o 

1.31467 

r * 

cH20 

55.4  mol  1_1 

ch2o 

-2.51444 

K 

4 

dh2o 

-1.72542 

°Me 

1.23  10-7  cm2  s'1 

7’crit 

1  Me 

512.5  K 

Pm? 

80.8  105  Pa 

-rent 

1  H2  O 

647.1  K 

Ph2o 

220.6  105  Pa 

Table  3 

Operation  conditions  used  for  simulation  of  all  scenarios. 


Anode: 

Cathode 

Parameter 

Value 

Parameter 

Value 

pa 

1.013  bar 

PC 

1.013  bar 

■jrA 

323.15  K 

f1 

323.15  K 

rha 

100% 

CMe.O 

1  mol  1_1 

cc 

o2,0 

7.92  10'3 

mol  1  1 

CA 

lh2o,o 

53.15  mol  1_1 

rc 

n2,0 

29.8  10'3 

mol  r1 

CA 

LC°2,0 

0  mol  r1 

CC 

co2,0 

0  mol  r1 

CA 

cMeg,0 

0  mol  1_1 

CC 

H2Og,0 

0  mol  1_1 

c a 

lh2os,o 

0  mol  r1 

Tc 

in 

293  K 

rA 

Me.in 

CA 

LH20,in 

1  mol  r1 

53.15  mol  1_1 

1.013  bar 
10 

rA 

LC02,in 

0  mol  1_1 

cc 

lC02  ,in 

0  mol  1_1 

CA 

LMeg,in 

0  mol  r1 

CA 

LH2Og,in 

0  mol  r1 
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scenario.  Furthermore,  it  is  assumed  that  the  gas  entering  the 
cathode  is  saturated  with  water.  Considering  these  assumptions, 
the  mass  balance,  Eq.  (8),  for  cathodic  water  in  steady  state  can  be 
reduced  to: 


0  =  ri 


H2Og,in 


rC  c  c  I 
cH2Ogrout  '  o 


H2Og 


(22) 


With  k  =  0  and  D^j1  =  0,  Eqs.  (19)— (21)  can  be  modified  as 
follows: 


AMi 
4 T 


(23) 


Aui 

If 


(24) 


(25) 


Including  Eq.  (12)  for  n^Qg  =  0,  Eqs.  (23)— (25),  (A.21),  (A.11)  and 

(A.12)  into  Eq.  (22)  results  in  an  equation  to  calculate  the  air  excess 
ratio  Ac  required  to  supply  sufficient  water.  A  detailed  derivation  of 
this  equation  can  be  found  in  the  appendix  (A.25): 


Ac 


2-3 


'HoOS 


'gas 


y 


H2Os,in 


'H-iOS 


'gas 


0.21  1 


(26) 


with  the  mole  fraction  of  water  at  the  cathodic  inlet  y^0  in  and  the 
water  concentration  Qg  at  the  cathodic  outlet.  For  ideal  gases,  the 
corresponding  norm  volume  flow  rate  at  cathodic  inlet  can  be 
calculated  as: 


RT] 


P  N 


nm 


(27) 


Including  Eqs.  (A.17)  and  (A.19)  results  in: 


RTn  1 

PN  0.21  (l-y^o, in) 


(28) 


Fig.  3.  Required  air  excess  ratio  for  sufficient  water  supply  by  cathodic  inlet  gas  as  a 
function  of  inlet  temperature  for  total  water  consumption  and  considering  diffusion 
through  a  GDL  of  thickness  dGDL  =  300  pm  at  a  current  density  of  i  =  200  mA  cm  2  and 
i  =  400  mA  cm  2  respectively. 


the  cathodic  gas.  Further  limitations  due  to  mass  transport  and 
dehumidification  kinetics  which  inhibit  total  dehumidification 
intensify  the  requirements.  This  is  indicated  by  the  grey  curves  in 
Fig.  3  which  include  water  vapour  transport  through  a  GDL  of 
300  pm  thickness  at  a  current  density  of  i  =  200  mA  cm-2  and 
i  =  400  mA  cm-2  respectively.  The  remaining  water  concentration 
in  the  cathode  chamber  c^0  is  then  calculated  by  Eq.  (A.22).  The 
results  show  that  high  inlet  temperatures  are  required  to  supply 
sufficient  water  for  reasonable  air  excess  ratios  if  mass  transport 
through  membrane  does  not  occur.  This  is  intensified  by  consid¬ 
ering  limitation  due  to  mass  transport  from  chamber  through  GDL 
to  reaction  region.  This  mass  transport  is  proportional  to  current 
density.  Thus,  if  the  transport  of  water  vapour  through  GDL  is 
considered,  higher  current  density  necessitates  higher  water 
diffusion  and,  thus,  requires  higher  water  vapour  concentration  at 
cathode.  Therefore,  higher  current  densities  lead  to  higher  air 
excess  ratios  required  for  sufficient  water  supply. 


Additional  parameters  needed  to  simulate  the  model  of  this 
scenario  are  listed  in  Table  4. 

4.1.2.  Results  and  discussion 

Since  this  scenario  is  only  considering  the  cathode  of  the  fuel 
cell,  the  results  shown  below  can  be  applied  to  all  alkaline  fuel  cells 
fed  with  air  at  the  cathode. 

Fig.  3  shows  the  minimum  air  excess  ratio  required  to  supply 
sufficient  water  at  the  cathode  as  a  function  of  inlet  temperature  for 
RHfn  =  100%.  The  limiting  curve  that  assumes  total  water  con¬ 
sumption  (c£20  =  0)  decreases  with  increasing  inlet  temperature. 
This  curve  is  independent  of  current  density.  Below  the  curve, 
water  supply  is  insufficient,  above  the  curve  some  water  remains  in 


Table  4 

Additional  parameters  used  for  simulation  of  scenario  1. 


Parameter 

Value 

Parameter 

Value 

dGDL 

300  pm 

i 

400  mA  cur2 

ogs 

0.15  cm2  s  1 

RHFn 

100% 

Tn 

273.15  K 

Pn 

1.013  bar 

Fig.  4.  Required  air  excess  ratio  for  sufficient  water  supply  by  cathodic  inlet  gas  as  a 
function  of  inlet  temperature  for  total  water  consumption  at  reference  conditions:  all 
gas  pressures  are  p  =  1  bar  and  RH^  =  100%;  and  at  p  =  2bar  as  well  as  for  RH^  =  60%. 
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Fig.  5.  Required  inlet  volume  flow  rate  of  air  for  sufficient  water  supply  by  cathodic 
inlet  gas  as  a  function  of  current  for  total  water  consumption  at  reference  conditions: 

=  300  K,  RH^  =  100%  and  p  =  1.013  bar,  and  at  7T  =  273  K,  at  p  =  2.026  bar  and  at 
RHpn  =  60%. 


Fig.  4  reveals  the  influence  of  pressure  and  relative  humidity.  All 
three  curves  in  this  figure  assume  total  water  consumption  at 
cathode.  The  reference  conditions  include  that  pressure  at  cathode 
and  cathodic  inlet  is  pc  =  pfn  =  p  =  1.013  bar  and  the  relative  hu¬ 
midity  at  cathodic  inlet  is  RH^n  =  100%.  These  conditions  are  in 
accordance  with  the  conditions  of  the  black  curve  in  Fig.  3.  A 
pressure  increase,  e.g.  top  =  2.026  bar,  necessitates  higher  values  of 
air  excess  ratio  in  order  to  supply  sufficient  water.  The  same  impact 
is  observed  if  the  relative  humidity  of  air  at  cathodic  inlet  is  low¬ 
ered,  e.g.  to  RHfn  =  60%.  That  is  because  in  both  cases  the  amount  of 
water  in  the  inlet  air  is  lower  and,  thus,  more  air  needs  to  be  fed  to 
the  fuel  cell  to  provide  sufficient  water. 

The  required  air  excess  ratios  in  Fig.  4  are  not  influenced  by 
current  produced  by  the  fuel  cell.  However,  the  corresponding  vol¬ 
ume  flow  rates  which  are  calculated  by  Eq.  (28)  for  the  same  con¬ 
ditions  do  increase  with  current  as  shown  in  Fig.  5.  A  norm  volume 
flow  rate  of  approximately  Vin  N  =  1.6/N  min1  is  required  to  supply 
sufficient  water  to  produce  a  current  of  /  =  10  A  at  reference  con¬ 
ditions.  A  decrease  of  inlet  temperature  to  =  273  I<  results  in  a 
volume  flow  rate  of  Vin  N  =  3/N  min'1 .  The  volume  flows  required  in 
the  cases  of  higher  pressure  (p  =  2.026  bar)  and  lower  humidity 
(RH^n  =  60%)  are  also  approximately  double  of  that  at  reference 
conditions.  Hence,  for  low  inlet  temperatures,  unreasonably  high 
volume  flow  rates  at  cathodic  inlet  are  required  for  high  currents, 
while  cathodic  water  supply  seems  feasible  for  low  currents. 

In  conclusion,  the  air  needs  to  be  humidified  and  either  high  inlet 
temperatures  or  high  volume  flow  rates  at  cathodic  inlet  are  required 
in  order  to  supply  sufficient  water  to  cathode.  These  two  options 
require  additional  components,  e.g.  humidifier  and  heater,  which 
significantly  complicate  the  fuel  cell  system  and  which  are  detri¬ 
mental  to  system  efficiency.  Hence,  this  scenario  reveals  that  water 
supply  by  cathodic  inlet  is  insufficient  and  mass  transport  through 
membrane  is  essential  for  cathodic  water  supply  in  alkaline  fuel  cells. 


Therefore,  the  general  model  is  slightly  modified  and  afterwards 
the  results  are  discussed. 


4.2.1.  Mathematical  modelling 

Water  is  supplied  at  cathode  by  oxidation  of  methanol  diffusing 
through  membrane  as  well  as  directly  by  water  diffusion  from 
anode  to  cathode.  In  order  to  analyse  the  requirements  on  water 
transport  through  membrane  for  sufficient  water  supply,  the  cor¬ 
responding  diffusion  coefficient  that  exactly  satisfies  cathodic 
water  demand  can  be  derived  from  Eq.  (8)  for  water  in  steady 
state: 


c  Fc 
H2Orout 


^H20,in 


dM 


AMcA 


(29) 


Instead  of  humidifying  cathodic  inlet  gas,  ambient  air  is  used  at 

=  293  K  and  RHfn  =  60%.  As  such,  gas  in  the  cathode  chamber 
does  not  supply  water  for  the  electrochemical  reaction.  On  the 
contrary,  it  is  assumed  that  either  some  of  the  available  water  is 
evaporated  and,  thus,  lost  for  reaction  or  the  water  content  of  the 
gas  does  not  change  at  all.  In  order  to  analyse  the  effect  of  electro- 
osmotic  water  drag,  methanol  cross-over  and  evaporation  at 
cathode,  these  processes  are  switched  off  one  by  one  and  the  re¬ 
sults  are  compared  to  the  reference  conditions  in  which  all  pro¬ 
cesses  occur.  Hence,  scenario  2  is  analysed  for  the  following 
conditions:  With  electro-osmotic  water  drag  k  =  4  or  without  k  =  0, 
with  methanol  cross-over  Dj^[e  =  1.23-10-7  cm2  s-1  or  without 
DMe  =  0,  with  evaporation  at  cathode  to  RHC  =  100%  at  7°  =  323  K, 
nH2o,out  >  "H2o,m  or  without  evaporation  ri^0  out  =  nfoin  and  for 
different  air  excess  ratios  Ac. 

Several  sources  and  sinks  of  water  change  the  amount  of  water 
at  anode.  In  the  following,  water  accumulation  AUH2o  denominates 
the  difference  between  the  volume  of  liquid  water  leaving  the 
anode  and  the  volume  of  liquid  water  fed  to  the  anode. 


AV 


h2o 


cA  cA 

rH20,out  rH20,in 


^H2Q,out  ~  ^H2Q,in 
CH20 


(30) 


Additional  parameters  needed  to  simulate  the  model  of  this 
scenario  are  listed  in  Table  5. 


4.2.2.  Results  and  discussion 

Fig.  6  displays  the  required  diffusion  coefficient  of  water  calcu¬ 
lated  by  Eq.  (29)  for  reference  conditions  in  which  water  drag, 
methanol  cross-over  and  evaporation  occur  and  Ac  =  10.  It  further 
shows  four  other  cases  which  differ  from  the  reference  condition  in 
one  of  the  above  mentioned  conditions  each. 

In  all  cases,  the  required  water  diffusion  coefficient  increases 
with  current  density  since  a  higher  current  leads  to  increased  water 
consumption.  Thus,  regarding  water  transport  through  membrane, 
the  cathodic  water  supply  is  mainly  a  challenge  for  high  current 
densities.  A  significant  influence  of  electro-osmotic  water  drag  and 
evaporation  can  be  observed  while  the  diffusion  coefficient  is 
insignificantly  higher  without  methanol  cross-over  for  a  methanol 
concentration  of  1  mol  1_1  at  anodic  inlet.  The  case  without  evap¬ 
oration  is  per  se  independent  of  air  excess  ratio.  In  the  case  that 


4.2.  Scenario  2  -  water  supply  by  mass  transport  through 
membrane 

This  scenario  is  used  to  study  the  requirements  on  the  water 
diffusion  through  membrane  for  sufficient  cathodic  water  supply. 


Table  5 

Additional  parameters  used  for  simulation  of  scenario  2. 


Parameter 

Value 

Parameter 

Value 

Aa 

4 

RHFn 

60% 

Ac 

10 

RHC 

100% 

C.  Weinzierl,  U.  Krewer  /  Journal  of  Power  Sources  268  (2014)  911—921 


917 


c 

0 

'o 


0 

o 

o 

c 

o 

CO 


"O 

L- 

CD 

CO 


§ 


A  A 

x  10~6 

\  .4 

—  Reference  conditions 

1.2 

.  Without  evaporation 

Without  water  drag  ^  * 

1 

■  Without  cross-over  S*  +  ^  ■ 

-  -  Oxygen  excess  ratio  Ac  =  4  S  s 

0.8 

/; ■ 

0.6 

S'' 

0.4 

S' 

0.2 

n 

^ _ 

0  ^  ■  ■  ■ 

0  100  200  300  400 

Current  density  [mA/cm2] 


Fig.  6.  Required  water  diffusion  coefficient  in  scenario  2  at  reference  conditions: 
Including  electro-osmotic  water  drag,  methanol  cross-over  and  evaporation  at  cathode 
with  an  air  excess  ratio  of  Ac  =  10;  and  for  four  other  cases  that  differ  from  reference 
conditions  in  one  of  the  mentioned  conditions. 


Current  density  [mA/cm2] 

Fig.  7.  Methanol  concentration  in  liquid  phase  at  anode  in  scenario  2  at  reference 
conditions:  Including  electro-osmotic  water  drag,  methanol  cross-over  and  evapora¬ 
tion  at  cathode  with  an  air  excess  ratio  of  Xc  =  10;  and  for  three  other  cases  that  differ 
from  reference  conditions  in  one  of  the  mentioned  conditions. 


evaporation  occurs,  the  required  water  diffusion  coefficient  de¬ 
creases  with  decreasing  air  excess  ratio.  Thus,  the  influence  of 
evaporation  is  reduced  in  case  of  a  smaller  air  excess  ratio. 

Mass  transport  through  GDL  and  CL  which  is  disregarded  in  this 
study  would  of  course  influence  these  results.  At  cathode,  a  diffu¬ 
sive  transport  of  water  from  CL  to  cathode  chamber  would  lead  to 
lower  evaporation  rates  and,  thus,  to  lower  required  diffusion  co¬ 
efficient.  Hence,  the  required  diffusion  coefficient  in  real  fuel  cells 
lies  in  between  the  cases  of  no  evaporation  and  full  evaporation.  At 
anode,  water  is  in  excess  and  is  produced  by  reaction  in  the  CL  but 
concentration  of  methanol  would  decrease  and  cross-over  would 
decrease  accordingly.  However,  the  influence  of  methanol  cross¬ 
over  on  the  required  water  diffusion  coefficient  through  mem¬ 
brane  is  insignificant. 

The  value  of  the  diffusion  coefficient  for  real  materials  is  no 
function  of  current  density.  Therefore,  the  simulated  values  of 
diffusion  coefficient  give  an  information  on  the  minimum  diffusion 
coefficient  a  membrane  needs  in  order  to  allow  operation  up  to  the 
given  current.  In  the  investigated  range  of  current  density,  the 
required  water  diffusion  coefficients  for  all  conditions  are  for 
example  far  below  those  observed  for  the  proton  exchange  mem¬ 
brane  Nation,  Dfj^(323  K)  =  4.6-10-6  cm2  s-1  [31.  The  only 
research  paper  about  water  diffusion  through  an  anion  exchange 
membrane  published  so  far  [12]  shows  diffusion  coefficients  close 
to  those  of  Nation  for  303  I<  and  313  K.  Hence,  it  may  be  expected 
that  the  water  diffusion  coefficient  through  an  anion  exchange 
membrane  is  sufficient  to  supply  sufficient  water  for  cathode 
reaction. 

Nevertheless,  there  will  be  a  discrepancy  between  the  required 
and  the  actual  diffusion  coefficient  of  water.  If  the  actual  diffusion 
coefficient  is  higher  than  the  required  one,  water  is  accumulated 
at  the  cathode  and  the  assumption  c^0  =  0  is  no  longer  valid. 
Consequently,  the  concentration  gradient  decreases  and  water 
diffusion  is  diminished.  Nevertheless,  it  is  possible  that  some 
pores  of  catalyst  layer  or  gas  diffusion  layer  at  the  cathode  get 
blocked  by  liquid  water  which  would  lead  to  a  decrease  in  cell 
performance. 

If  the  actual  diffusion  coefficient  is  below  the  required  one, 
water  transport  through  membrane  is  not  sufficient  for  cathodic 
water  supply.  Either  the  lacking  water  has  to  be  taken  from  the  gas 


or  the  current  density  has  to  be  reduced  until  the  water  diffusion  is 
sufficient. 

The  resulting  methanol  concentrations  in  the  liquid  phase  at 
anode  as  a  function  of  current  density  are  displayed  in  Fig.  7.  In  the 
case  without  methanol  cross-over,  methanol  concentration  does 
not  change  with  current  density.  This  can  be  explained  as  follows: 
All  inlet  volume  flow  rates  are  proportional  to  current  density  due 
to  their  definition  by  constant  excess  ratios  (see  Eqs.  (A.20)  and 
(A.21)).  Sources  and  sink  terms  consisting  of  reaction  rates  and 
the  electro-osmotic  drag  of  water  (Eqs.  (13)— (21))  are  also  pro¬ 
portional  to  current  density.  Since  water  diffusion  is  defined  to  be 
equal  to  water  consumption  at  cathode,  water  diffusion  is  also 
proportional  to  current  density.  Hence,  all  terms  needed  to  calcu¬ 
late  outlet  flow  rates  in  Eqs.  (11)  and  (12)  are  proportional  to  cur¬ 
rent  density  as  well.  Consequently,  all  terms  of  the  molar  balances 
(Eqs.  (7)  and  (8))  in  steady  state  with  rip  in  =  Cp  inFin  and 
rijg  out  =  CpF0Ut  are  proportional  to  current  density  which  can  be 


Fig.  8.  Methanol  efficiency  in  scenario  2  with  a  methanol  excess  ratio  of  XA  =  4  at 
reference  conditions  with  liquid  recycling,  without  cross-over  and  without  liquid 
recycling.  Reference  conditions:  Including  electro-osmotic  water  drag,  methanol  cross¬ 
over  and  evaporation  at  cathode  with  an  air  excess  ratio  of  Xc  =  10. 
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cancelled.  This  results  in  a  system  of  algebraic  equations  that  is 
independent  of  current  density  and  can  be  solved  for  steady  state 
concentrations  of  components.  Thus,  all  steady  state  concentrations 
are  a  function  of  excess  ratios,  inlet  concentrations,  concentrations 
of  pure  substances,  water  drag  coefficient,  pressure  and  tempera¬ 
ture  at  anode  and  cathode  but  are  independent  of  current  density. 
However,  the  concentration  of  methanol  in  liquid  phase  at  anode  is 
below  the  inlet  concentration,  since  methanol  is  consumed  while 
water  is  produced  by  electrochemical  reaction.  Since  methanol 
concentration  is  independent  of  current  density,  methanol  effi¬ 
ciency  is  also  independent  of  current  density  and  is  nearly  100%  as 
shown  in  Fig.  8.  The  slight  efficiency  loss  is  caused  by  methanol 
evaporation  at  anode. 

If  methanol  cross-over  occurs,  sources  and  sinks  at  cathode  and, 
consequently,  water  diffusion  and  outlet  volume  flows  contain 
terms  that  are  not  proportional  to  current  density  and  concentra¬ 
tions  are  no  longer  independent  of  current  density.  The  higher  the 
current  density,  the  smaller  is  the  influence  of  the  term  indepen¬ 
dent  of  current  density.  As  a  consequence,  the  influence  of  meth¬ 
anol  loss  by  cross-over  is  substantial  at  low  current  densities  and 
causes  small  methanol  concentrations  and  low  efficiencies.  At  high 
current  densities,  methanol  loss  caused  by  cross-over  is  small 
compared  to  the  other  processes  and  methanol  concentration 
converges  to  the  value  of  the  case  without  cross-over  for  high 
current  densities. 

In  the  three  cases  that  include  cross-over,  methanol  concen¬ 
tration  differs  slightly  from  that  at  reference  conditions  since  the 
required  water  diffusion  is  smaller  in  case  of  lower  air  excess  ratio 
Ac  =  4  or  without  evaporation  and  methanol  dilution  is  more 
intense.  However,  the  difference  is  small  and  the  corresponding 
change  in  efficiency  is  insignificant.  Hence,  the  methanol  effi¬ 
ciency  in  Fig.  8  is  only  analysed  for  three  different  cases.  The  first 
two  cases  show  methanol  efficiency  in  reference  conditions  and 
without  methanol  cross-over  assuming  that  all  methanol  at  anodic 
exit  can  be  recycled.  The  third  case  shows  methanol  efficiency 
without  methanol  recycling.  Methanol  cross-over  causes  high  ef¬ 
ficiency  losses,  especially  for  low  current  densities.  Furthermore, 
high  fuel  efficiencies  can  only  be  achieved  by  recycling  the 
methanol  at  anodic  outlet.  Since  the  methanol  concentration  for 
all  cases  shown  in  Fig.  7  is  below  the  inlet  concentration,  recycling 


Fig.  9.  Accumulation  of  liquid  water  at  anode  in  scenario  2  at  reference  conditions: 
Including  electro-osmotic  water  drag,  methanol  cross-over  and  evaporation  at  cathode 
with  an  air  excess  ratio  of  Ac  =  10;  and  for  three  other  cases  that  differ  from  reference 
conditions  in  one  of  the  mentioned  conditions. 


of  the  methanol  solution  requires  treatment  to  keep  the  methanol 
concentration  constant  and  to  avoid  dilution. 

The  amount  of  water  in  the  recycled  liquid  also  needs  to  be 
considered.  Fig.  9  shows  the  water  accumulation  at  the  anode 
which  is  calculated  by  Eq.  (30).  Both,  negative  and  positive  water 
accumulation  are  possible.  In  the  case  without  evaporation,  the 
required  water  diffusion  is  quite  small  and,  therefore,  more  water  is 
produced  at  anode  than  removed  by  diffusion.  As  a  result,  the  water 
accumulation  is  positive  because  the  amount  of  water  leaving  the 
anode  is  bigger  than  the  amount  of  water  fed  to  the  anode.  All  other 
cases  lead  to  a  negative  water  accumulation  which  means  that  the 
molar  flow  of  water  at  the  anodic  inlet  is  bigger  than  at  the  outlet. 
In  these  cases,  the  fuel  cell  runs  out  of  water.  This  indicates  that 
recycling  the  liquid  phase  at  the  anode  will  lead  to  another 
important  challenge  regarding  water  management  in  ADMFCs:  the 
water  level  stabilisation. 

5.  Conclusions 

As  in  most  kinds  of  fuel  cells,  water  management  is  a  chal¬ 
lenging  topic  in  alkaline  direct  methanol  fuel  cells.  Water  is 
consumed  at  the  gas-flushed  cathode  and  produced  at  the  anode 
where  water  is  already  in  excess.  In  this  paper,  two  scenarios  of  an 
ADMFC  with  different  operation  conditions  are  modelled  and 
analysed  in  order  to  study  conditions  for  sufficient  water  supply  to 
cathode.  Water  management  demands  as  well  as  the  influence  of 
processes  like  water  drag  and  methanol  cross-over  on  water 
management  in  ADMFCs  are  identified. 

Scenario  1  reveals  that  humidifying  cathodic  inlet  gas  is  not 
sufficient  for  water  supply  to  cathode  and  it  is  important  to  consider 
water  management  at  anode  and  water  transport  through  mem¬ 
brane  for  any  alkaline  anion  exchange  membrane  fuel  cell.  If  water  is 
only  supplied  by  cathodic  inlet,  the  air  fed  to  cathode  needs  to  be 
saturated  with  water  and  either  high  inlet  temperatures  or  high 
volume  flow  rates  are  required  to  supply  sufficient  water.  This  is 
manageable  for  small  systems  operated  in  the  lab  for  analysis,  but  it 
is  unrealisable  for  real  fuel  cell  systems  that  are  supposed  to  work 
with  high  efficiency.  Hence,  the  simulations  reveal  that  water 
transport  through  membrane  is  essential  and  desired  to  supply 
sufficient  water  to  the  cathode  for  all  alkaline  anion  exchange 
membrane  fuel  cells. 

Scenario  2  shows  that  the  required  diffusion  coefficient  of 
water  for  sufficient  water  supply  to  the  cathode  of  an  ADMFC 
depends  on  operation  conditions  and  rises  with  increasing  current 
density.  For  membranes  with  a  diffusion  coefficient  of  water  close 
to  that  of  Nation,  water  diffusion  is  sufficient  to  achieve  current 
densities  of  i  <  400  mA  cm-2.  Thus,  additional  humidification  of 
cathodic  inlet  gas  is  necessary  at  high  current  densities,  but  for 
low  current  densities,  humidification  might  cause  flooding  of 
cathode.  Therefore,  lack  of  water  is  not  the  reason  for  the  low 
performance  of  present  ADMFCs  but  flooding  might  be.  The  results 
also  show  the  necessity  to  consider  the  diffusion  coefficient  of 
water  when  designing  new  membranes  for  alkaline  fuel  cells  since 
water  diffusivity  is  fixed  by  membrane  material.  The  influence  of 
methanol  cross-over  on  the  requirements  placed  on  the  mem¬ 
brane  are  negligible  but  electro-osmotic  water  drag  and  water 
evaporation  at  cathode  have  a  significant  effect.  The  latter  depends 
on  relative  humidity  of  cathodic  inlet  gas,  air  excess  ratio  and  the 
equilibrium  of  humidity  of  gas  phase  with  the  moistness  of  the 
membrane. 

Stable  water  management  at  cathode  of  an  ADMFC  can  be 
achieved  by  adjusting  evaporation  and  current  density. 
Furthermore,  the  operation  range  can  be  influenced  by  air 
excess  ratio,  water  content  of  inlet  air  and  also  by  membrane 
material. 
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In  order  to  achieve  high  methanol  efficiencies,  it  is  necessary  to 
recycle  methanol  solution  at  the  anodic  outlet.  However,  it  is  not 
sufficient  to  consider  methanol  content  in  the  recycled  liquid  but 
also  water  content  is  important  since  the  fuel  cell  can  run  out  of 
water  or  accumulate  water  depending  on  the  operation  conditions. 
Thus,  water  level  stabilisation  at  the  anode  is  another  challenging 
topic  for  further  studies. 


The  concentrations  of  liquid  methanol  c^e  and  water  0  at 
anode  are  based  on  the  volume  of  the  whole  anode  chamber  which 
includes  gas  and  liquid  phase.  For  calculating  diffusion  of  methanol 
and  water  the  concentrations  based  on  liquid  volume  only  are 
needed.  These  can  be  calculated  as  follows: 


Me,liq 
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Me 


VA 

iiq 


cA  VA 

LMev 

VA 

iiq 
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Appendix  A.  Additional  model  equations 

The  following  additional  equations  are  used  for  the  modelling  of 
all  scenarios  described  in  this  paper: 
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The  flow  rates  entering  or  leaving  anode  or  cathode  are 
described  by: 
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In  Scenario  1  (Section  4.1)  the  remaining  concentration  of  water 
in  the  cathode  chamber  is  calculated  by  including  a  diffusion  layer. 
The  diffusive  molar  flow  of  water  is  equal  to  the  water  consumption 
by  reaction.  Therefore,  the  minimal  required  water  concentration 
can  be  calculated  by: 
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Thereby,  the  vapour  pressure  in  Eq.  (A.9)  is  described  by  the 
Wagner  equation  [30]. 
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This  concentration  is  also  needed  for  the  calculation  of  the 
required  air  excess  ratio  Ac  in  Eq.  (26)  that  is  derived  in  the 
following: 


Including  Eqs.  (23)— (25)  in  Eq.  (12)  for  riJ^Qg  =  0  results  in: 
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Combining  Eqs.  (A.21 ),  (A.ll )  and  (A.12),  the  molar  flux  of  water 
entering  the  cathode  chamber  is  calculated  by: 
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For  ideal  gases,  concentrations  of  pure  gases  are  equal  at  the 
same  conditions  =  c^Qg  =  c^as  and  Eq.  (26)  can  be  derived  by 
combining  Eqs.  (A.23),  (A.24)  and  (22): 


Nomenclature 


Latin  symbols 

Am  area  of  membrane  [m2] 

A  Wagner  coefficient  for  vapour  pressure  [-] 

B  Wagner  coefficient  for  vapour  pressure  [-] 

C  Wagner  coefficient  for  vapour  pressure  [-] 

Cp  concentration  of  component  f3  [mol  l-1] 

concentration  of  pure  substance  (3  [mol  1_1] 

D  Wagner  coefficient  for  vapour  pressure  [-] 

D¥  diffusion  coefficient  of  component  (3  through  membrane 
[m2  s'1] 

dM  thickness  of  membrane  [m] 

F  Faraday  constant  [A  s  mor1] 

Fin/out  volume  flow  rate  at  inlet/outlet  [m3  s'1] 

i  current  density  [A  m'2] 

jk,p  diffusion  flux  of  (3  in  direction  k  [kg  m'2  s] 

lip  molar  flow  of  component  (3  [mol  s'1] 

p  pressure  [Pa] 

pcnt  critical  pressure  [Pa] 

p°(F)  vapour  pressure  of  (3  at  temperature  T  [Pa] 

RH  relative  humidity  [%] 

T  temperature  [K] 

Tcnt  critical  temperature  [K] 

Tr  temperature  referred  to  critical  temperature  [-] 

t  time  [s] 

tmax  maximal  simulated  operation  time  [s] 

V  volume  [m3] 

v/<  velocity  in  direction  k  [m  s'1] 

y  mole  fraction  in  gas  phase  [-] 


X  mole  fraction  in  liquid  phase  [-] 

Zk  space  coordinate  in  direction  k  [m] 

Greek  symbols 

k  water  drag  coefficient  [-] 

A  excess  ratio  of  reactant  [-] 

pp  density  of  component  (3  [mol  l'1] 

(Tp  mass  sources  and  sinks  of  component  / 3  [kg  s'1] 

cr p  molar  sources  and  sinks  of  component  j3  [mol  s'1] 

Superscripts 
A  at  anode 

C  at  cathode 

co  cross-over 

crit  critical  value 

diff  diffusion 

drag  electro-osmotic  water  drag 
exc  excess 

M  of  membrane 

Sys  system 

Subscripts 

a  component  in  anode  chamber 

j3  component  general 

7  component  in  cathode  chamber 

dry  dry  gas 

k  spatial  direction 

N  at  norm  conditions 

in  at  inlet 

liq  of  liquid  phase 

out  at  outlet 

Me  liquid  methanol 

Meg  methanol  vapour 

C02  carbon  dioxide 

H20  liquid  water 

H2Og  water  vapour 

N2  nitrogen 
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O2  oxygen 


Abbreviations 

AAEM  alkaline  anion  exchange  membrane 

ADMFC  alkaline  direct  methanol  fuel  cell 

AFC  alkaline  fuel  cell 

CSTR  continuous  stirred  tank  reactor 

FC  fuel  cell 


References 

[1]  E.H.  Yu,  U.  Krewer,  K.  Scott,  Energies  3  (2010)  1499-1528. 

[2]  E.H.  Yu,  X.  Wang,  U.  Krewer,  L.  Li,  K.  Scott,  Energy  Environ.  Sci.  5  (2012) 
5668-5680. 

[3]  A.  Tewari,  V.  Sambhy,  M.  Macdonald,  A.  Sen,  J.  Power  Sources  153  (2006) 
1-10. 

[4]  H.  Ko,  H.  Juang,  J.  Appl.  Electrochem.  13  (1983)  725—730. 

[5]  J.  Varcoe,  R.  Slade,  E.  Lam  How  Yee,  Chem.  Commun.  (2006)  1428—1429. 

[6]  L.A.  Adams,  S.D.  Poynton,  C.  Tamain,  R.C.T.  Slade,  J.R.  Varcoe,  ChemSusChem  1 
(2008)  79-81. 

[7]  E.  Yu,  K.  Scott,  J.  Power  Sources  137  (2004)  248-256. 

[8]  L.  Sun,  J.  Guo,  J.  Zhou,  Q.  Xu,  D.  Chu,  R.  Chen,  J.  Power  Sources  202  (2012) 
70-77. 

[9]  N.  Fujiwara,  Z.  Siroma,  S.-i.  Yamazaki,  T.  Ioroi,  H.  Senoh,  K.  Yasuda,  J.  Power 
Sources  185  (2008)  621-626. 

[10]  J.R.  Varcoe,  M.  Beillard,  D.M.  Halepoto,  J.P.  Kizewski,  S.D.  Poynton,  R.C.T.  Slade, 
ECS  Trans.  16  (2008)  1819-1834. 


[11]  S.D.  Poynton,  J.P.  Kizewski,  R.C.T.  Slade,  J.R.  Varcoe,  Solid  State  Ionics  181 
(2010)  219—222,  14th  International  Conference  on  Solid  State  Protonic  Con¬ 
ductors  (SSPC-14),  Kyoto,  Japan,  Sep  07—11,  2008. 

[12]  Y.S.  Li,  T.S.  Zhao,  W.W.  Yang,  Int.  J.  Hydrogen  Energy  35  (2010)  5656-5665. 

[13]  N.  Follain,  S.  Roualdes,  S.  Marais,  J.  Frugier,  M.  Reinholdt,  J.  Durand,  J.  Phys. 
Chem.  C  116  (2012)  8510-8522. 

[14]  T.D.  Myles,  A.M.  Kiss,  K.N.  Grew,  A.A.  Peracchio,  G.J.  Nelson,  W.K.S.  Chiu, 
J.  Electrochem.  Soc.  158  (2011)  B790-B796. 

[15]  T.  Yamanaka,  T.  Takeguchi,  H.  Takahashi,  W.  Ueda,  J.  Electrochem.  Soc.  156 
(2009)  B831-B835. 

[16]  I.  Verhaert,  S.  Verhelst,  G.  Janssen,  G.  Mulder,  M.  De  Paepe,  Int.  J.  Hydrogen 
Energy  36  (2011)  11011-11024. 

[17]  S.  Rowshanzamir,  M.  Kazemeini,  M.  Isfahani,  Int.  J.  Hydrogen  Energy  23  (1998) 
499-506. 

[18]  P.  Bjornbom,  S.  Yang,  Electrochem.  Acta  38  (1993)  2599—2609. 

[19]  S.  Yang,  P.  Bjornbom,  Electrochem.  Acta  37  (1992)  1831—1843. 

[20]  S.  Huo,  H.  Deng,  Y.  Chang,  K.  Jiao,  Int.  J.  Hydrogen  Energy  37  (2012) 
18389-18402. 

[21]  H.  Deng,  S.  Huo,  Y.  Chang,  Y.  Zhou,  K.  Jiao,  Int.  J.  Hydrogen  Energy  38  (2013) 
6509-6525. 

[22]  M.  Kimble,  R.  White,  J.  Electrochem.  Soc.  138  (1991)  3370—3382. 

[23]  J.  Jo,  S.  Moon,  S.  Yi,  J.  Appl.  Electrochem.  30  (2000)  1023—1031. 

[24]  S.  Mohan,  S.O.B.  Shrestha,  J.  Fuel  Cell  Sci.  Technol.  7  (2010). 

[25]  I.  Verhaert,  M.  De  Paepe,  G.  Mulder,  J.  Power  Sources  193  (2009)  233—240. 

[26]  A.  Verma,  S.  Basu,  J.  Power  Sources  168  (2007)  200—210. 

[27]  K.N.  Grew,  W.K.S.  Chiu,  J.  Electrochem.  Soc.  157  (2010)  B327-B337. 

[28]  K.N.  Grew,  D.  Chu,  W.K.S.  Chiu,  J.  Electrochem.  Soc.  157  (2010)  B1024-B1032. 

[29]  C.Y.  Du,  T.S.  Zhao,  W.W.  Yang,  Electrochem.  Acta  52  (2007)  5266-5271. 

[30]  Verein  Deutscher  Ingenieure,  VDI-Gesellschaft  Verfahrenstechnik  und  Chemie 
Ingenieurwesen  GVC  (Ed.),  VDI-Warmeatlas,  Springer-Verlag,  Berlin,  Heidel¬ 
berg,  New  York,  2006. 

[31]  S.  Yeo,  A.  Eisenberg,  J.  Appl.  Polym.  Sci.  21  (1977)  875—898. 


